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Abstract 

We discuss the bremsstrahlung of photons into a heat bath, and calculate from first principles the 
energy radiated. Even to lowest order the spectrum of the radiation at low frequency is no more 
singular than at zero temperature. In addition to the obvious contributions, this spectrum includes 
■ terms associated with fluctuations. 



u . 

In a recent paper, Weldon W has considered the spectrum of photons emitted when a charged particle 
\ scatters off a potential within a heat bath. As he observes, in such a problem infrared effects are more 

delicate than at zero temperature, because in addition to the usual factors the integrands involved 
, contain Bose distributions 

£ ! ^ = m 

. which diverge at zero energy. According to Weldon, while the infrared divergences do nevertheless 
still cancel in a lowest-order perturbative calculation (this is an extension of the zero-temperature 
Bloch-Nordsieck mechanism), they leave behind a spectrum of emitted photons whose number density 

. dN/du) at zero-frequency u> is unacceptable: it behaves as u>~ 2 , so that the energy density udN/du 
behaves as which integrates to give a divergent result for the total energy carried off by the 

photons*. 

Q_i| Weldon further shows how to sum the infrared-sensitive terms to all orders. In this note, we are not 
concerned with this, the most important part of Weldon's beautiful paper. Rather, we seek to clarify 
his remarks about the lowest-order calculation of the radiated energy. We do so within the framework 
.£h ! °f a very simple mo dell 2 ! which we believe incorporates the essential features. This is the decay, 
within a heat bath that contains only photons, of a scalar uncharged particle P into a pair of scalar 
charged particles ("muons"). We calculate three different quantities: (I) the order-e 2 contribution 
Ar to the total decay rate; (II) the probability that the muons emerge with combined centre-of-mass 
energy within AE of the parent-particle energy; and (III) the net energy spectrum ojdN/duj transferred 
to the heat bath by the photons. 

For (I), we agree with Weldon that Ar is finite. For (II), we find that the probability in order e 2 
behaves as (AE) -1 ; we expect that Weldon's resummation procedure would remove this singular 
behaviour. It is for (III) that we differ from Weldon: ujdN/duo is finite at u = even to order e 2 . 

As is well-known, the decay rate may be calculated from a discontinuity of a thermal Green's 
function PI. However, it appears that the energy spectrum udN/duj cannot be calulated in terms 
of such a discontinuity. Instead, we must introduce a different method, which we derive from first 
principles. This method reproduces the same result for the decay rate as the conventional one, and 
may be used also for calculating oodN/doj. 

We recall first that in finite-temperature gauge theory one may choose that at bare-propagator level 
only the physical degrees of freedom are thermalisedW . Then the thermal part of the photon propa- 
gator 

A T (k) = 2ir5(k 2 )f(k) (2a) 



* See equation (1.7) of reference 1 



represents the effects of real photons in the heat bath. This is easy to understand. The finite- 
temperature photon propagator is, in any gauge, 



Z~hi e- Ha/T TA(x)A(0) = Z~hr e~ Ha/T ({Q\TA(x)A(0)\Q)+ : A(x)A(0) : ) (3a) 

where 

Z = tr e~ Ho/T (36) 

and the traces are defined as summation over expectation values in only physical states. The first 
term in (3a) is just the zero-temperature Feynman propagator in whatever gauge is being used, and is 
unaltered by the heat bath. The second receives contributions only from the real photons in the heat 
bath and is 

Z~hi e~ H ^ T [ -^-27T5 + (k 2 )-^-2Tr5 + (k' 2 ) fat (k)a(k')e ik - x + J \k')a{k)e- ik A (26) 

J (27TJ 4 (27TJ 4 I ) 

which is the Fourier transform of (2a). 

We return now to the decay in the heat bath of the uncharged particle into two muons. We have 
previously^] calculated, to lowest order e 2 in the muon charge, the change Ar in the decay rate 
caused by the heat bath. This change involves only the thermal part (2) of the photon propagator. 
It consists of three pieces: (i) a term that may be interpreted as the acquisition by the muons of an 
additional mass of order eT, (ii) another that may be interpreted as a shift induced by the heat bath 
in the coupling A that governs the decay, and (iii) a term that corresponds to stimulated emission of 
photons into the heat bath or absorption from it. We found that the mass shift (i) is 

5m 2 = 2e 2 J -0^2nS(k 2 )f(k°) (4a) 

and the terms (ii) and (iii) contribute 

2P° Ar = AV^L. J # h 5(k 2 )f(k^ Pl d*p 2 5 + (p 2 - m 2 )5 + {p 2 - m 2 ) - ( J?L - JEL) 

{-5\pi + p 2 - P) + 5\pi + P2 + k - P)} (46) 

where in the last factor the first 5-function corresponds to the term (ii) and the second to (iii). We 
have supposed that the decaying particle has momentum P. 

In the integrals (2), k denotes photon 4- momentum. Because only the thermal part (2) of the photon 
propagator is involved, the photon is real, even when its effect is interpreted as just changing the value 
of the coupling A, or of the muon mass. These two changes together correspond to the absorption of 
a real photon from the heat bath and its re-emission with the same momentum. The amplitude for 
this is the sum of the order-e 2 graphs shown in figure la and lb, where the dashed lines at the bottom 
are the spectator photons in the heat bath. Both of the active photons in the graphs must have the 
same momentum k, because to obtain an order-e 2 contribution to the decay rate we must multiply 
the sum of graphs in figures la and lb by the order-e graph of figure lc, in which all the heat-bath 
photons are spectators. This product corresponds to the first two of the of cut graphs shown in figure 
2, in which the photon line represents the thermal part (2a) of the photon propagator. Because this 
thermal propagator contains 5(k 2 ), rather than 5 + (k 2 ), each single cut graph of figure 2 corresponds 
to two contributions, as indicated in in figures la and lb. The first <5-function in the curly bracket in 
(4b) arises from the first two cut graphs in figure 2, while the second (5-function arise from the other 
two cut graphs. The latter represent the sum of the squares of the terms shown in figure Id and le. 

In his calculation of the net photon energy spectrum todN/duj, Weldon appears not to have included 
the contributions of figure la, lb and le. We maintain that it is necessary to do so if one wishes to 
calculate a physically-meaningful quantity. In order to construct such a quantity, we must carefully 
perform a thought experiment. We shall consider two such experiments. In both, we think of the heat 
bath as a microwave oven with only photons in it. Any muons that are produced as a result of decays 



escape straight through the oven walls. Our world has no other charged particles, in particular no 
electrons. 

In our first thought experiment, we drill a small hole in one of the oven walls and imagine observing 
the energy spectrum codN/duj of radiation emitted through it, just as one does in the usual analysis 
of black-body radiation^. We suppose that we have introduced into the oven a collection of the 
particles P that decay into two muons, and consider the spectrum of photons in events triggered by 
the decay-product muons. We calculate the difference between this photon spectrum at temperature 
T and that at zero temperature. 

In order to do this calculation, we go back to first principles. Define single-photon states r, initially 
discrete with energies e r , and label the states of the heat bath by the corresponding occupation 
numbers n r . Then, with a decaying particle P also in the oven, the initial density matrix is 

Pi = Z' 1 ^exp (- ^n r e r /T) \n 1 ,n 2 , . . . ,P)(n 1 ,n 2 , . . . ,P\ (5a) 
W r 

with 

r 

After a long time, the density matrix has become 

Pf = Z' 1 ^exp (- ^n r e r /T) T\m,n 2 , ■ ■ .,P)(n 1 ,n 2 , . . -,P\T^ (6a) 
{"} r 

Here, as usual, T is the interaction part of the 5-matrix, S = 1 + iT. Because our initial state contains 
an unstable particle, which inevitably decays, we know that an interaction takes place. The decay 
rate of the particle P is calculated from the expectation value of pf in the states 

\m 1 ,m 2 ,...[x(p 1 )[x(p 2 )) (66) 

with summation over the photon occupation numbers mi,m 2 , . . . and over the momenta p\ and p 2 
of the decay-product muons. For this, we need the photon-operator structure of the T-matrix. The 
various matrix elements in figure 1 correspond to the following photon-operator terms in T: 

e 2 A a a\a a + e 2 ^ B s a s al + C + e ^ D s a\ + e ^ E s a s (7) 

s s s s 

The operators A S ,B S , . . . correspond to figures la,b,. . . and involve the photon modes s only through 
conserving energy and momentum. They convert the initial particle P into the two muons and 
include also the internal muon propagators. We use the familiar properties a\n) = y/n \n — 1) and 

at|n) = y/n + l \n + 1), so that the order-e 2 contribution to the summed expectation value takes the 
form 

Ar = Z- 1 ^ exp ( - ^ n r e r /T) ^ (n s P s AC + (n s + 1)P S BC + (n s + 1)P S DD + n s P s EE ) (8a) 

{n} r s, Pl ,p 2 

Here 

Ps AC = e 2 (p(pi)p( P2 )\ A s \P)(P\ Ct \p( Pl )p(p 2 )) + e 2 (p( Pl )p( P2 )\ C \P)(P\ a\ \p( Pl )p(p 2 )) (86) 

and similarly for P EC \P® D and Pf E . Notice that, in this order-e 2 expression for Ar, the DD and 
EE terms are the squares of order-e amplitudes, while the AC and BC terms are the interference 
between order-e 2 and order-e amplitudes. The presence of the interference terms at finite temperature 
is possible because the order-e amplitude in figure lc involves spectator photons in the heat bath. 
One of these spectator photons has to be identified with the the photon that is absorbed and emitted 
in the diagrams of figure la and lb. For this to be possible, the absorbed and emitted photon must 



(56) 



have the same momentum as the corresponding spectator photon in figure lc, and so the photon 
scattering in figures la and lb is forward scattering. 

We use 

n = Z- 1 ^ne-" e / T = /(e) (9) 

n 

where /(e) is the Bose distribution (1), so that (8a) is 

Ar= (f(e s )P* C + (f(ts) + l)P s BC + (f(e s ) + l)P? D + f(e s )P s EE ) (10) 

S,Pl,P2 

Subtracting off the zero-temperature limit, we retrieve in the continuum limit the formula (4b). 

We now calculate the thermal average of the number of photons in state t. This comes from the 

expectation value of a\a t pf in the states (6b), again with summation over the photon occupation 
numbers mi, m 2 , ■ ■ ■ and over the momenta p\ and p 2 of the decay-product muons. Because pf is not 
normalised to unity, we must introduce also a factor 



(N t ) = T-'Z- 1 Y, exp ( - Y n rtr/T) Y { n t P t AC + Mnt + l)P t BC + K + l) 2 P t DD + (n t - l)n t P ( 

{n} r Pi,P2 



EE 
t 



+n t P cc + n t Y (n s P s A ° + K + l)^f ° + K + 1)P? D + n.P? E ) } (11) 



We use 

so that (11) is 
where 



n 2 



= Z~ X Y n 2 e~ ne/T = n 2 + /(e) (/(e) + l) (12) 



t 



?EE 

t 



(N t ) = f(e t ) + F emission (e t ) + F absorption (et) + F fluc (e t ) (13a) 

F emission (et)=r -l ^ (/( et ) + l)p/ 

Pl,P2 

F absorption (e t ) = -r- 1 Y f( e t)Pi 

Pl,P2 

F fluc (e t )=r- 1 Y f(et){f{tt) + l)[Pf C + P t BC +P t DD + P t EE ] (136) 

Pl,P2 

Evidently, the first term in (13a) is the value (Nt) would have if there was no decaying particle in the 
heat bath. We subtract it off. As expected, the emission term, corresponding to part of the square 
of figure Id, adds to the spectrum, while the absorption term, corresponding to part of the square 
of figure le, has negative sign and so reduces it. From (12), the term F fluc (e t ), because it involves 
/(/ + 1), is associated with fluctuations in the heat bath. 

Subtracting off the zero-temperature limit, and going over to the continuum limit, we are left with 
the change in the spectrum per decaying particle: 

2P° T—AN(u) = A emission + A absorption + A fluc (14a) 

with 

A 2 p 2 r 

^mission + ^bsorption = /( w )^^ / d 4 k5(k° - Co)5(k 2 )d 4 Pl d 4 p 2 5 + (v\ ~ m 2 )5 + ( V \ ~ m 2 ) 

(2vr)° J 



J 2 {^{Pi +p 2 + k-P)- 5\ Pl +p 2 -k-P)} 



(146) 



and 



A fluc = /M(/M + J d ^k5(k° - cv)S(k 2 )d 4 Pl d 4 P 2S + (p 2 - m 2 )5 + {p 2 - m 2 ) 

J 2 {s 4 (Pi +P2 + k-P) + 6 4 ( Pl + P2 -k-P)- 25\ Pl +P2-P)} 

A\ 2 p 2 f 8 1" 

+/M(/M + 1)— J d 4 k5(k° - ^(k 2 )^ J d 4 Pl d 4 p 2 5 + (p 2 - m 2 )5 + (p 2 - m 2 )5\ Pl + p 2 -P) 

(14c) 

where 

j2 = _ ( ll_ _ P2 \ 

\P!.k p 2 .k) 

The last term in (14c) is the part of the interference between figures la and lb with figure lc that is 
associated with the thermal mass shift 5m 2 of (4a). 

When lj — ► 0, -£^AN(lo) diverges as lo -1 . The uj~ 2 mentioned by WeldonW is present in i? absor P tlon 
and F emisslon ^ but it cancels between them and leaves behind a; -1 , The separate contributions to the 
fluctuation term behave as u>~ 3 . An important component of the fluctuation term is the interference 
between the diagrams of figures la and lb with figure lc; the Sm 2 term is a part of this interference 
ternJ 2 ]. It is straightforward to see that the co~ 3 behaviour cancels between the various parts of the 
fluctuation term, as does its oj~ 2 behaviour. Moreover, one can show that, no matter how fast the 
decaying particle is moving (or not moving) before it decays, the w _1 behaviour cancels also; this 
cancellation is closely related to the cancellation that occurs at low temperature in Ar, making^ 
it behave as T 4 instead of T 2 . As we show in the Appendix, it is a result of current conservation 
together with unitarity. So the fluctuation term is nonsingular when to — ► 0; it actually goes to zero 
like oj. However, notice that the uj~ x behaviour in the fluctuation term cancels only because we have 
integrated over the momenta p\ and p 2 of the final-state charged particles. In the scattering problem 
considered by WeldonM, pi rather is the momentum of an incoming particle and is not integrated, so 
we believe that then the cancellation of the w _1 contributions to the fluctuation terms does not take 
place. 

Thus the breakdown of perturbation theory signalled by Weldon in fact does not occur for the mea- 
surement he considered. However, it does occur in certain other cases. To see this, consider our second 
thought experiment (II), where instead of measuring the number density of the photons emitted from 
the heat bath we measure the total energy of the two muons. We consider the rate at which the pair 
of muons emerges with total energy within magnitude AE of the energy of their parent particle. This 
is calculated from (4b) with an additional 

e(AE-\P Q -p\-pl\) (15) 

under the integral. In spite of the presence in (4b) of the factor that is quadratically divergent at 
k = 0, the fc-integration is convergent, because although the difference of (5-functions in the last factor 
is linear in k, the rest of the integrand is even in k. However, although this cancellation of divergences 
takes place between the two 5-functions, as AE — ► the answer is again divergent. The term with 
the first (5-function clearly is independent of AE, while for small AE the second behaves as AE -1 . 
This behaviour is unphysical (it is even large negative) and so must be corrected by a resummation 
such as was suggested by Weldon W. Of course, the same resummation is necessary for this quantity 
at zero temperature, though then before the resummation the divergence is only as log AE. 

We have explained that our method of calculation requires us to go back to first principles. It gives 
the same answer for the decay rate as the conventional approach, but may be applied also to the 
calculation of the photon energy spectrum ujdN/doj or the number spectrum dN/duj. While the decay 
rate is linear in the Bose distribution f(u), the number spectrum dN/dco contains a quadratic term, 
which is associated with fluctuations in the heat bath. Therefore, it does not seem to be associated 
with a discontinuity of a thermal Green's function. Our initial attempts to calculate dN/duj from the 



(Ud) 



discontinuity diagrams of figure 2 yielded an absorption term whose sign was opposite from that in 
(13b), which is clearly incorrect; as we have shown, the correct sign is obtained only when one takes 
account of the fact that part of the absorption effect is a contribution to the fluctuation term. 
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Appendix 

In this Appendix, we calculate the leading terms in the photon spectrum (14) for small u. We verify 
that the oo~ 2 terms cancel in the emission and absorption parts of (14b), leaving u~ l . We find further 
that there is a cancellation between the different parts of the fluctuation term (14c), so that while 
the separate terms behave as u~ 3 their sum behaves as uj. This is closely related to the cancellation, 
found in [2], of the low-temperature contributions to AI\ As we explain at the end of this Appendix, 
it is a consequence of current conservation and unitarity. 

In order to deal with (14), we follow reference 2 and make the changes of variables 

Pi-k Pi-k 

Pl ~* Pl T "pT ' P2 ^ P2Zf ~pj k ( A1 ) 

so that 

J 2 5 + (p\-m 2 )5 + (p 2 2 -m 2 )5 A {p 1 +p 2 -P±k) — ► J 2 5 + {p\ - m 2 T X)5 + {p\ - m 2 T X)5 A ( Pl +p 2 - P) 

(A2) 

where 

X = 2pi t^ k (A3) 
P.k 

We then expand (A2) in powers of X, to the second order for the fluctuation contribution to (14), 
and to first order for the emission and absorption contributions. The result is 

\ 2 p 2 r 

A emission + ^bsorption _ f{u,)±l- / d A k5{k 2 )5(k° - Uj)A (A4a) 

(2vr) & J 

and 

A fluc ~ /M(/M + 1)^ J d 4 k5(k 2 )5(k° - u)B (A4b) 
A = j ' d i Pl d i p 2 5 i {P - Pl -p 2 )(-2XJ 2 )(d/dm 2 )5 + (p 2 -m 2 )5 + (p 2 -m 2 ) (A5a) 

B = j d A Pl d 4 p 2 5 4 (P - Pl -p 2 ) [X 2 J 2 (d/dm 2 ) 2 + 4(d/dm 2 )]5 + (pl-m 2 )5 + (p 2 2 -m 2 ), (A5b) 
with J 2 defined in (14d). 

Evidently both A and B are Lorentz invariants and so they can depend upon k via the invariant P.k 
only. Therefore we can first calculate with the P-particle at rest (relative to the heat bath), and then 
replace Mlo by P.k to get the result for a moving P-particle. For P at rest, and with 6 the angle 
between k and pi, 

1 f 1 u vt2 2M T 2m 2 , (M + J(M 2 -4m 2 )\l /A , 
- / d(cos0) AJ 2 = 1 — — In ^ — ;jf A6a) 

2 7_i ' uj [ M^(M 2 - Am 2 ) \M - ^(M 2 - Am 2 ) )\ V ; 



where 



and 



- f 1 dxX 2 J 2 = -(M 2 -4m 2 
2 J -i 3 



{A6b) 



We have anticipated that we are going to insert these integrals into (A5) and (A4), and so set k 2 = 
and p\ = m 2 =p\. We require also 

J d* Pl d*p 2 6*(P- Pl -p 2 )5+(p 2 - m 2 )5+(p 2 - m 2 ) = ff - 4m 2 ) ^ 

Since (A6b) and (A6c) do not depend on u, they remain the same when the decaying particle P is not 
at rest in the heat bath. Using (A6) we find that the two derivatives in the expression (A5b) for B 
give equal and opposite contributions when they are inserted into (A4b). This cancellation is exactly 
similar to that which occurs in the change Ar induced by the heat bath in the low-temperature decay 
ratet 2 !. 

Thus A fluc vanishes in this approximation: it receives its first contribution from the X 4 terms in the 
expansion of (A2) and so behaves as w 3 f(oj)(f(uj) + 1) at low frequency. The leading behaviour of 
—i 



-^AA f (w) near uj = comes from the emission and absorption terms and is, from (A5a) and (A6a), 



2p o T dAN 2X 2 e 2 /H / M + ^(M 2 - 4m 2 ) \ / P° + V(P° 2 - M 2 ) \ 

duo (2vr) 3 ^(po2 _ M 2) \M - J{M 2 -Am 2 )) \p° - J (P* 2 - M*) ) K ' 

where we have used 



/ 



d 4 kS(k° - ^{^{P.k)- 1 = In f P0 + ^ P0 ] M 2 " 



^(P° z - M 2 ) \P°- - M 2 ) 



The reason for the cancellation in Ar, and of the leading terms in B, is the following. Ar is calculated 
from the trace of the final-state density matrix pf given in (6a). Through unitarity, this is 

tr Z~ x ex P ( — 5^ n r e r /T) T\ni,n 2 , ■ ■ ■ , p)(m,n 2 ,...,P|rt 

{»} 

= -iZ~ l exp (- n r e r /T) (m, n 2 , . . . , P\(T - T^)\m, n 2 , . . . , P) (A8) 
{»} r 

Because we are working to lowest order in the charge e, only one photon in the state |ni, n 2 , . . . , P) is 
active, with the other photons just providing the number factors n s in (8a), so that in the continuum 
limit Ar is an integral, involving the Bose factor /, over the photon momentum k of the imaginary 
part of the scattering amplitude (P, k^\T\P, k v ). Because the electromagnetic current is conserved, 
for off-shell k this amplitude has the usual decomposition: 

-(^ - k -Th + ( pti - ii k ")( p " - ^ 

This has to be contracted with transverse real-photon polarisation vectors and k 2 set to 0. This 
eliminates T 2 . But when k 2 — ► 0, Xi ~ (P.k) 2 T 2 /k 2 , and T 2 vanishes like k 2 . Because we do not 
allow the unstable particle P to survive in the final state, the amplitudes T\ and T 2 do not have 
poles corresponding to this particle, and so T\ behaves as (PA;) 2 for small P.k. It is this that makes 
Ar behave^] as T 4 at low temperature T. According to (8a), when the zero-temperature limit has 
been subtracted off, Ar is calculated from the combination [P t AC + P t BC + P t DD + P t EE ] . This same 
combination appears in the fluctuation term (13b) and the (P.k) 2 there gives the additional factor uj 2 . 



Note that it is essential to this argument that we have integrated over the momenta p\ and p 2 of the 
final-state muons. 
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Figure 1 

Photon radiation in a heat bath, involving both active and spectator photons. 




Figure 2 



Cut diagrams for the square of the amplitude corresponding to the various terms of figure 1 



